GLOBAL WELL-POSEDNESS FOR CUBIC NLS WITH 
NONLINEAR DAMPING 
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Abstract. We study the Cauchy problem for the cubic nonlinear Schrodinger 
equation, perturbed by (higher order) dissipative nonlinearities. We prove 
global in-time existence of solutions for general initial data in the energy space. 
In particular wc treat the energy-critical case of a quintic dissipation in three 
space dimensions. 



1. Introduction 

In this paper we study the Cauchy problem of the cubic nonlinear Schrodinger 
equation (NLS) with nonlinear damping, i.e. 

J id t u+^-Au=V(x)u + X\u\ 2 u-ia\u\ p ~ 1 u, (t, x) € [0, oo) x R d , 
(1.1) < 2 

[ u(0) =u (x), 

with given parameters A € K and a ^ 0, the latter describing the strength of the 
dissipation within our model. We shall consider the physically relevant situation 
of d ^ 3 spatial dimensions and assume that the dissipative nonlinearity is at least 
of the same order as the cubic one, i.e. p 3. However, in dimension d — 3, we 
shall restrict ourselves to 3 ^ p ^ 5. In other words, we assume that the dissipative 
effects in our model can be described by nonlinearities, which are at most energy 
critical in the sense of [12, 35]. In the following, we shall also assume the external 
potential V to be an anisotropic quadratic confinement, i.e. 



(1.2) V(x) = ±J2^xl Uj G M. 

j=i 

In the case without any dissipation, i.e. a = 0, equation (1.1) simplifies to the 
classical cubic NLS, a canonical description for weakly nonlinear wave propagation 
in dispersive media [34]. It arises in various areas of physics, such as nonlinear 
optics, plasma physics or ultracold quantum gases (Bose-Einstein condensates), cf. 
[34] for a broader introduction. Due to the inclusion of a quadratic potential, the 
natural energy space when studying the Cauchy problem for this equation is given 
by 

S:={ii£ i7 1 (M rf ) : xu £ L 2 (R d )} , 
and we consequently denote 

INe := IMU* + HVuHls + \\xu\\ L 2. 
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In the case without harmonic confinement, the results on local and global well- 
posedness properties of (purely dispersive) cubic NLS are by now considered to be 
classical, see e.g. [12, 35] and the references given therein. In particular, it is very 
well known that in the focusing case A < finite-time blow-up of solutions can 
occur if d > 2, i.e. 

lim ||Vu(i,-)||z,2( R d) = oo, 

where the blow-up time T < oo depends on the initial data uq. Generalization 
to the case with harmonic potentials have been studied in [17, 26, 38] and more 
recently in [7, 9, 10, 11]. For confining potentials, the well-posedness results are 
found to be very much the same as in the classical situation (see also the discussion 
in Remark 4.3 below). 

From the point of view of physics, the occurrence of blow-up usually implies that 
new effects have to be taken into account in order to extend the model beyond the 
appearing singularities. As a specific example, let us briefly describe such a situ- 
ation in the context of collapsing Bose-Einstein condensates: These are ultracold 
quantum gases, which (in a mean-field approximation) can be described via the 
Gross-Pitaevskii equation, a cubic NLS governing the macroscopic wave function of 
the condensate. In addition, one usually takes into account a harmonic confinement 
V, modeling the electromagnetic trap present in actual physical experiments. It is 
nowadays possible to study the collapse of the condensate experimentally by tuning 
the effective nonlinear interaction from positive to negative. What one observes is, 
that as the particle density increases around the blow-up point, atoms are sud- 
denly emitted from the condensate in bursts, so-called jets. These jets are caused 
by (inelastic) three-body interaction or recombination effects taking place only at 
high densities [1] . Mathematically, three-body forces can be effectively described by 
quintic nonlinearities, cf. [13]. In order to describe the emittance of particles from 
the condensate within the realm of Gross-Pitaevskii theory one usually employs the 
following dissipative model (see [1, 21, 31] and the references given therein): 

h 2 

(1.3) ihd t u + - — Au — V{x)u + g\u\ 2 u — ia^ufu. 

Here g = 4irNH 2 a/m, with m denoting the mass of the particles and a their scat- 
tering length. The parameter 03 denotes the three particle recombination loss-rate. 
Rescaling (1.3) into dimensionless form yields (1.1) with p = 5 and d = 3. Note 
that in (1.3) the total mass, i.e. M{t) := \\u(t, OIII2, is no longer conserved, as can 
easily be seen from the dissipation equation for the particle density p = \u\ 2 : 

(1.4) d t p + Aw J= -2ap 3 , 

where J := Im(uVu) denotes the current density. Numerical simulations show 
[5] (see also [4]), that the solution to (1.3) undergoes a series of events where 
concentration is followed by consequent dissipation, lowering the density and hence 
preventing blow-up. We consequently expect that equation (1.3) and more generally 
(1.1), admits global in-time solutions (even for large initial u € S) and it is the 
main purpose of our work to prove that this is indeed the case. Obviously, we are 
mainly interested in the focusing case A < 0, but we shall keep A e K for the sake of 
generality. Let us also remark that, even though equation (1.3) has been the main 
motivation for the present work, similar NLS type models with nonlinear damping 
terms also appear in other areas of physics, sec e.g. [6, 28, 29, 32]. For example, 
in the context of nonlinear optics, nonlinear damping terms are used to describe 
multi-photon absorption, cf. [16]. 

From a mathematical point of view, we are facing two basic problems in setting 
up a global existence theory for (1.1): First, in the physically most relevant situa- 
tion where p = 5 and d = 3, the (quintic) dissipative nonlinearity is known to be 
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energy-critical [12, 35] and thus cannot be considered as a small perturbation of 
the Laplacian any more. Indeed, global well-posedness for defocusing energy crit- 
ical NLS in d = 3 dimensions (without external potentials) has only been proved 
recently in the seminal work of Colliander et al. [14] (see also [23]). Even though 
our nonlinearity is dissipative, this property is not seen locally in time (e.g. in 
terms of Strichartz-estimates) and we are therefore in a similar situation as for the 
Hamiltonian case. A second obstruction for a global existence-theory, even for en- 
ergy sub-critical cases, is the lack of conserved quantities, in contrast to the usual 
case of Hamiltonian NLS. This is also the main difference to the work of Tao et al. 
[36], which treats purely dispersive NLS with combined power type nonlinearitics. 
Also there, blow-up is prevented by sufficiently strong (i.e. higher order) defocus- 
ing nonlinearities. Note however, that such a (mass conservative) model would not 
be able to describe the coherent loss of particles found in collapsing Bosc-Einstcin 
condensates. 

Before going further, let us briefly compare our situation to the mathematically 
much better studied case of linearly damped NLS, i.e. 

(1.5) id t u + ^Am = V(x)u + X\u\ 2 u - iau. 

There exists a considerable amount of results in the physics and mathematics lit- 
erature for such (weakly damped) NLS type equations. In particular the Cauchy 
problem to (1.5) has been analyzed in [27, 37] and more detailed properties con- 
cerning the long-time behavior of solutions can be found in e.g. [16, 19, 25]. The 
main difference between our case and (1.5) is, that the latter can be treated by the 
phase-transformation u(t,x) — ► e at u(t,x), which makes the damping term vanish. 

In contrast to the situation with linear damping, the literature on NLS with 
nonlinear damping is not so abundant: In [24, 33] the asymptotic behavior of small 
solutions to NLS with dissipative nonlinearitics of the form A|u| p_1 u, with Im A 
and 1 < p ^ 1 + 2/d, is studied. For higher order dissipative nonlinearitics, the 
only rigorous results, we are aware of, are given in [28], where, for p = 2(1 + s) with 
s > 0, the authors prove a non-uniform (in-time) bound on the H 1 norm of the 
solution. Numerical studies of NLS type models with nonlinear damping can be 
found in [4, 5] and we also mention the results of [16], based on modulation theory. 

In the following section we shall present our main theorem and deduce from it 
several corollaries. The corresponding proofs are then given in Section 3 (where the 
energy critical case is treated) and Section 4. 



Our first main result concerns the case of an energy-critical damping term. 

Theorem 2.1. Let d = 3 and p = 5. Assume V to be a quadratic confinement of 
the type (1.2) and m € E. Then, for any A e M and for any a > the equation 
(1.1) has a unique global in-time solution u G C([0, oo),E), such that 



The theorem confirms the numerical results of [5] on the time-evolution of at- 
tractive Bose-Einstein condensates with three-body recombination. No collapse, 
i.e. finite-time blow-up, occurs due to the dissipation. Note that Theorem 2.1 
holds for any a > 0, no matter how small, which is consistent with the physics 
literature, where, after a dimensionless rescaling, one finds a <C 1 for any realistic 
parameter regime [1, 5, 21]. The additional a-priori estimate on the Ll° x norm of 
the solution (to be proved in Proposition 3.2) is reminiscent of the one used in [14]. 
This space-time bound is also required in our proof, since it is well known (sec e.g. 



2. Main results 
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[12]) that the usual a-priori bound on the H 1 norm is not be sufficient to conclude 
global existence. The reason being that the local existence-time of solutions does 
not only depend on the H 1 norm of u, but also on its profile. 

As a first consequence of the above given theorem, we obtain the analogous state- 
ment for all energy subcritical situations (in d < 3) where the nonlinear damping 
term is of higher order than cubic. 

Corollary 2.2. Let V be a quadratic confinement and iio £ S. Moreover, let 
p > 3, if d = 1, 2 and 3 < p < 5, if d = 3. Then, for any A € R and any a > 0, the 
equation (1.1) has a unique global in-time solution u e C([0,oo),£). 

The result strengthens the one of [28] , as we are able to prove uniform (in-time) 
a-priori bounds on the energy. As a by-product of our analysis we also obtain a 
well-posendess result in the case where the damping term is not of higher order but 
given by a cubic nonlincarity, i.e. we consider 

(2.1) id t u + ^Au = V(x)u + (A - ia)\u\ 2 u. 

In the context of Bose-Einstein condensates, equation (2.1) corresponds to a model 
where one only takes into account two-body losses, see [1, 31]. Most of the time, 
however, they are neglected in view of the three-body recombination effects. 

Corollary 2.3. Let d ^ 3, assume V to be quadratic and u € £. Then, for 
any a ^ max{0, — A} the equation (2.1) has a unique global in-time solution u £ 
C([0,oo),E). 

Thus, even in the case where the two nonlinearities exactly balance, i.e. a = |A|, 
we have a global in-time solution. This indicates that dissipation acts on a faster 
time-scale than the focusing nonlincarity does. In situations where A < (focussing 
case) and a < |A| we expect that solutions to (2.1) in general exhibit finite-time 
blow-up (see also the discussion in Remark 4.3). It remains an interesting open 
problem to rigorously prove that this is indeed the case. 

Remark 2.4. In the case where V(x) = 0, equation (2.1) is similar to the complex 
Ginzburg- Landau equation, which usually reads 

(2.2) d t u = (1 + ia)Au - (1 + i(i)\u\ 2 u + 

for some given parameters a, /3, 7 G R, see e.g. [2] for a broader introduction. 
In the case of no external driving field 7 = 0, the main difference between the 
nonlincarly damped cubic NLS (2.1) and the Ginzburg-Landau equation (2.2) is 
that the latter invokes an additional linear diffusion oc Au. Equation (2.1) can 
therefore be considered as the diffusionlcss limit of (2.2) and our well-poscdncss 
results can consequently be reinterpreted within this limiting regime of the complex 
Ginzburg-Landau equation. For an application of this model in the context of 
nonlinear optics see [32]. 

We finally state the following result on time-decay of solutions as t — > +00. 

Corollary 2.5. Let u e C([0,oo),£), be a global-in-time solution to (1.1). Then 
u(t) decays to zero as t — > +00, in the following sense: For each sequence of time- 
steps (t n )neN tending to infinity, there exists a subsequence (t ne ), such that 

(2.3) lim \\u(t ne )\\ L 2 (Rd) = 0. 

I— >oo 

Moreover, let u n (t,x) := u(t + t n ,x), then, up to extraction of subsequences, we 
also have u n -» in L 2 (0, T); ij£ c (R d )). 
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Unfortunately, we are not able so far to derive explicit rates for the time-decay 
of the (total mass of the) solution. Note that the decay rates obtained in [24, 33] in 
general do not hold in our case, since we do not restrict ourselves to small solutions. 

Remark 2.6. In the fluid dynamics picture of quantum mechanics the equation for 
the mass density p — \u\ 2 has to be complemented by an equation for J = Im(uVu), 
forming the so-called quantum hydrodynamic system: 

( d t p + div J= -2ap( p+1 ^ 2 , 

< 2 - 4 > ( dtJ+ div (£2£) + , w + * v - i ,v (M) - w-^j. 

Observe that the equation for the current density J also picks up a nonlinear damp- 
ing term. The connection between (1.1) and (2.4) can be established following 
[3, 18], in order to translate our well-posedness results for (1.1) into the analogous 
statements for (2.4). 

3. The energy critical case 

This Section is devoted to the Proof of Theorem 2.1. That is, we want to show 
that solutions to 

(3.1) id t u+ ^Au = V(x)u + A|ii| 2 u- ia\u\ A u, (t, x) e [0, oo) x R 3 , 

with w(0) = u (x) e E, exist globally in-time. To this end, we shall first derive 
several a-priori bounds on the solution. 

3.1. A-priori estimates. From the dissipation equation (1.4), we immediately 
obtain the following lemma. 

Lemma 3.1. Let u(t) € E be a solution of (3.1). Then it holds 

(3.2) \\u(t,-)\\ L 2^\\u \\ L 2, Vt^O, 
and in addition we have 

poo 

(3.3) / \\u(t,-)\\ 6 L6 dt^C(\\u \\ L 2). 

Jo 

Proof. We argue formally and first integrate (1.4) w.r.t. x to obtain 

d f 
-^\Ht,-)\\h=-2o- / \u{t,x)\ 6 dx^0, 
at J M 3 

which consequently yields (3.2). If we then also integrate w.r.t. t, we get 
2a [ [ \u(t,x)\ 6 dxdt = [ \u(0, x)\ 2 dx - [ \u(t, x)\ 2 dx < f \u(0,x)\ 2 dx 

JO JR3 JR3 JR3 J R 3 

and thus (3.3). By a standard density argument this a-priori estimates are easily 
shown to hold for any u(t, •) € E. □ 

To proceed further consider the following energy-type functional 

(3.4) E K {t):= I l\Vu(t,x)\ 2 + V(x)\u(t,x)\ 2 + ^-\u(t,x)\ 4 + K \u{t,x)\ 6 dx, 

for some parameter k > to be chosen later on. The following a-priori bound on 
this functional will be a key ingredient in the proof of Theorem 2.1. 
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Proposition 3.2. Let u(t) € S be a solution of (3.1) and V(x) a quadratic con- 
finement of the form (1.2). Moreover, let < k < cr/6. T/ien, it ZioWs 

E K (t)^E K {0) + C(\\u \\ L 2), Vt>0, 

w/iere C = C(k, |A|, <t) > and we afeo Ziaue £/ie following space-time bounds 

f f \u(t,x)\ 4 \Vu{t,x)\ 2 dxdt ^C{E K (0),\\u \\ L 2), 

JO Jr 3 

f f V(x)\u(t,x)\ 6 dxdt^C(E K (0),\\u \\ L 2). 
Jo Jr 3 

In particular, it also holds 

[ [ \u(t,x)\ 10 dxdt^ C(E K (0),\\u o \\ L2 ). 
Jo Jr 3 

Proof. Let us assume u is regular enough to perform the following formal manip- 
ulations. A standard density argument can then be used to justify all of them 
rigorously. Computing the time-derivative of E K (t), we obtain 

d E K (t) = a[ M 4 Rc(uAu) dx - 2a [ V\u\ 6 dx - 2a\ [ \u\ s dx 
Jr 3 Jr 3 Jr 3 



dt 



— 3k |ii| 4 Im(uAw) dx — &Ka / |u| 10 cix 
Jr 3 Jr 3 

-a p\Vp\ 2 dx-a / p 2 |Vu| 2 dx - 2a\ \ p 4 dx 

Jr 3 Jr 3 Jr 3 



3k 



R3 



p Im(uAu) dx — 6(tk / p dx — 2a / V(x)p dx 



R3 



= -a [ p\\7p\ 2 dx-a [ p 2 \\7u\ 2 dx-2a\ [ p 4 dx 
Jr 3 Jr 3 Jr 3 

+ 6 k P^p ■ Im(wVii) dx — 6aK / p 5 dx — 2a 

Jr 3 Jr 3 Jr 3 



V(x)p 3 dx, 



where for the last equality we perfom an integration by parts, using the fact that 
Im(SAu) = div (Im(ltVu)). Next, we note that Vp = 2Re(uVu) and thus 



p|Vu| 2 = |uVw| 2 = |Rc(uVu)| 2 + |Im(uVu) 



J|Vp| 2 + |J| 2 , 



where, as before, J = Im(uVu) denotes the current density. Using this, we can 
rewrite 

2 



/ pVp • Im(i/Vti) dx = — p 
Jr 3 Jr 3 

Jr 3 



1. 



-Vp - J 



^Vp-J 



This consequently yields 



(3.5) 



a [ p|Vpp 

Jr 3 



dx 



dx 



•J 

Jr 3 



/pQivpi 2 +ui 2 

p 2 \Vu\ 2 dx. 



dx 



— 6k 



1. 



dx — (a — 6k) J p |Vu| dx — 2a\ j p dx 

2 



-Vp - J 



dx — 6aK / p 5 dx — a V(x)p 3 
Jr 3 Jr 3 



dx. 



Under the assumption k < a/6, all the terms on the right hand side are negative 
(recall that V(x) is assumed to be a quadratic confinement), except for the one 
proportional to ||p||l 4 ) which we shall treat by interpolation. Since 



\L 4 



,5/8 
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we have 



// dx < ^- / p 3 dx + ^- I p 5 dx, 
2e 7b3 2 



for some arbitrary constants e > 0. Choosing e < 6k/|A|, wc consequently obtain 
^-E K (t) < -cr /* p\Vp\ 2 dx-(a-6K) [ p 2 \Vu\ 2 dx-cr [ V(x)p 3 dx 



(3.6) 



■6k / P ^ v /°- J dx-Ci(e) [ p 5 dx + C 2 (e) [ p 3 dx. 

7r3 2 J R 3 J R 3 



There still remains a positive term on the r.h.s., namely the last one. However, we 
already know from (3.3) that 



/ / p 3 dxdt^C(\\u \\ L 2). 
Jo Jr 3 



We can therefore integrate (3.6) w.r.t. time and using the fact that E K (0) < oo by 
assumption, we obtain the assertion of the proposition, provided k < a/6. □ 

Remark 3.3. The proof shows, that it would not be sufficient to consider only the 
energy functional of the unperturbed equation, i.e. 

Eo(t)= I \\Vu{t,x)\ 2 + V(x)u{t,x) + ^-\u{t,x)\ A . 

In fact, computing the time-derivative of E (t), we find 

-^-E (t) = - [ p\Vp\ 2 dx- [ p 2 \Vu\ 2 dx-2 [ Vp 3 dx-2\[ p 4 dx. 

dt J R 3 J K 3 J R 3 J R 3 

Thus, if the A < (focusing cubic nonlinearity) the last term on the r.h.s. is positive 
and hence we can not conclude that Eo(t) is non increasing. 

As an immediate corollary of Proposition 3.2 we obtain the uniform boundedness 
of the energy norm of u(t). 

Corollary 3.4. Let V be quadratic and u(t) a solution to (3.1). Then 

||u(t,-)|| E <(7(||uo||j:), Vi)0. 

Proof. We already know from Lemma 3.1 that the L 2 norm of u(t) is bounded. If 
A > then the uniform bound on E K (t) immediately yields the assertion. On the 
other hand, if > A = — |A|, we write 

||Vu(t, + \\xu(t, .)||| 2 < 2E (t) + ^\\u(t, .)||l4, 
since V is quadratic. Now, by interpolation \\u\\\ 4 < IMIl 2 IMIisi and thus 

||V«(t, + \\xu{t, < 2E (t) + l -^\\u(t, -)\\ 2 L2 + ^||u(i, -)||!e 
for some e > 0. If we then choose e = grn, we get 

||V«(t, Oil! 2 + WMt, -)\\h < 2E K (t) + —\Ht, -)\\h 

Hi 

2X 2 

^ 2E K (0) + ||«o|| ia + C(||uo|UO- 

Since the left hand side is the sum of two non-negative term, each of them is 
bounded individually by a uniform constant C, depending only on ||u ||s- O 
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3.2. Strichartz estimates. In order to prove global well-posedness of (3.1) we 
shall heavily rely on the use of Strichartz estimates. Let us briefly recall the defini- 
tion and main properties of these estimates for the following Schrodinger propagator 

U(t) = e~ itH , H:=~A + V(x), 

where V is given by (1.2). The operator U(t) consequently generates the linear time- 
evolution corresponding to (1.1). We first note that in view of Mehler's formula, cf. 
[11, 7], the group U{t) is not only bounded on L 2 (R d ), but also enjoys dispersive 
properties for small time. More precisely it holds 

C 

(3-7) H^W/llz/^) < up7^ll/IUi(M<*) ; for 1*1 < 5, 

for some S = 6(ojj) > 0, see [7, 8]. Note that for an harmonic potential as in (1.2), 6 
is necessarily finite, since H has eigenvalues (see also below) . From this dispersive 
estimates one is led to the following (local in-time) Strichartz estimates for U(t) in 
terms of admissible index pairs. 

Definition 3.5. A pair (q, r) is admissible if 2 ^ r ^ (resp. 2 < r ^ oo, if 
d = 1 and 2 ^ r < oo, if d = 2) and 

2 -=d(l- 1 - 



q \2 r, 

Then, for any space-time slab / x R d , we can define the Strichartz norm 

H/lls°(7xR d ) : = SU P H/llz,«L£(IxR<')) 

(q>r) 

where the supremum is taken over all admissible pairs of exponents (q, r). 

Following [22] (see also [7, 8]) one has the following estimates, where (q',r') 
denotes the Holder dual exponents of (q, r): 

Lemma 3.6. Let (q,r), (qi,ri) and (q2,r 2 ) be admissible pairs. Let L be some 
finite time interval. Then it holds 

\\U(-M\ Lq{I:Lr) ^C(r,d)\I\^y\\ L2 , 

and also 



U(t-s)F(s)ds ^C(r u r 2 ,d)\I\^\\F\\ L9Hl;Lr ^ 



ln{s<t} 

Proof. This result can be essentially be found in [8, Propositon 3.3], but for the 
convenience of the reader we shall recall the basic idea in the case of an isotropic 
confinement, where the Hamiltonian is imply given by 

H = \^\*\ 

with oj G R. In this case, Mehler's formula yields local in-time dispersion on 
\t\ < 7r/(2w). Now, let J be a given finite time-interval. Then we can split I into 
(finitely many) sub-intervals Ij, j = 1, . .., N such that \Ij\ < n/(2w). Strichartz 
estimates (based on Mehler's formula) imply that on each Lj it holds 

\\u(tM Lq{LnLn ^cM L 2, 

for a universal C = C(r, d) > 0. From this we obtain 

N 
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This directly yields the first assertion of the lemma for any finite time-interval / in 
the case of an isotropic potential. A similar argument can then be done in order 
to prove the second assertion stated above and a generalization to the case of a 
non-isotropic confinement is straightforward. □ 

Taking ip to be an eigcnfunction of the (anisotropic) harmonic oscillator shows, 
that in general one can not expect the above given Strichartz-estimates to hold 
uniformly in-time, i.e. without dependence on the length of /, unless all the luj in 
(1.2) are in fact zero, in which case the operator U(t) simplifies to the usual free 
Schrodinger group. In the upcoming section, this requires us to keep track of the 
dependence of all appearing constants on the length of I. 

3.3. Proof of Theorem 2.1. Having in mind the a-priori bounds obtained in 
Section 3.1, we can now state the proof of our main result. 

Proof. We rewrite (3.1) using Duhamel's formula 

(3.8) u(t) = U{t)u + i\ [ U(t-s)(\u\ 2 u)(s)ds- a f U(t - s) (\u\ 4 u) (s) ds. 

Jo Jo 

We aim to prove global well-posedness by a fixed point argument. To this end, we 
first consider a space-time slab I x R 3 , where \I\ < +oo, such that the L\° x norm 
of u within this slab is small, say 

(3.9) IMIliojjxrs) < £ < 1- 

Then, using Strichartz and Holder estimates, we can estimate (3.8) as follows: 




and taking into account the smallncss assumption stated above we get 

H«IUf iS (/xRS) < (||«0|U= + |/| 1/2 £ 2 |M| S o + £ 4 |M| S o) • 

Next, in order to bound Vw we first note that 

[d j ,H]=d j V(x), [xj,H] = dj, j = l,...,d. 

where [A, B] = AB — BA denotes the usual commutator. By assumption, djV(x) = 
ujjXj, i.e. linear in x. This shows that we can obtain a closed family of estimates 
for Vu and xu. More precisely, we have 

l|Vu|| L | L r + |MU?££ ^ l 7 | 1/9 (l|Vu ||i,2 + \\xu n \\ L 2 + IHu^Vullii^ + ||a;|u| 2 u|| L i L 2 

+ III U | 4VU IIl2 l 6/5 + \\X\U\ 4 U\\ L 2 L 6/^ 

< l/lV^HVuollso + \\xuo\\ S o + IllVVfllVullso + \\xu\\ S o) 
+ e\\\\7u\\ s0 + \\xu\\ s0 )). 
Thus, denoting the Strichartz norm in S by 

Ms* ■= \\u\\so + \\Vu\\so + \\xu\\ S o, 

we infer 




Thus, if s < 1 defined in (3.9), is sufficiently small, a standard contraction argument 
yields 

Hs s <<7(||uo|| E ,|J|). 
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Next, consider any finite time interval I = [0, T], for some T < +00. From 
Proposition 3.2 we already know that the Lj° x of u is uniformly bounded but not 
necessarily small, say ||u||^io (jxrs) ^ M, where M is independent of the length of I. 
Then we can divide I into subintervals I = TiU. . .U/jv, where \h\ — C(S, ||u||e) > 0, 
and such that the L\° x norm is sufficiently small in each It = [tt-i,tt], i.e. 

IMIlJ^xr 3 ) < for alii = 1,. ..,N. 

Note that for any I = [0, T] this e < 1 only depends on ||uo||s and not on the length 
of the interval I, in view of Proposition 3.2. 

By the same fixed point argument as before, we consequently obtain that in each 
I t x R 3 it holds 

IMI&KJixR") < C(||«(t/_i, -)||e, |/«|), £ = 1, . . . ,N. 

By Corollary 3.4 we also have that \\u(t, is uniformly bounded for alH > and 
thus 

IMIs E (/ f xR 3 ) C(||« ||e, I^D- 
Summing up all the subintervals It we consequently infer 

IMIs s (/xr 3 ) «S C"(||u ||e,M). 

By continuity, we consequently obtain a unique solution u in [0, T] x R 3 , for any 
T e (0,oo) and thus we conclude that the Cauchy problem (3.1) is globally well- 
posed in S (depending continuously on the initial data). □ 

As a by-product of our analysis we infer that all Strichartz norms are uniformly 
bounded during the time evolution, at least in situations without harmonic con- 
finement. 

Proposition 3.7. Let d = 3, p = 5 and assume V(x) = 0. Then for any admissible 
pair of exponents (q, r) it holds: 

ll M lli9([0,oo):L'-(R 3 )) + l|Vu|| L g([0 : oo);L''(R 3 )) ^ C( 1 1 u 1 1 L 2 ! E K (0)) . 

Proof. Let (q, r) be an arbitrary admissible pair of exponents, I an arbitrary time 
interval (which could also be infinite), and let t* G I. Then, by Strichartz estimates 
we have 

l|Vu|| L fz, S (i X R3) < ||Vu(r)|| L 2 (R 3 ) + |||u| 2 Vu|| L ioA (/xR3) + |||u| 4 Vu|| i io/7 (/xR3) , 

where we recall that (10/3, 10/3) is an admissible pair of exponents. By Holder's 
inequality we obtain 

IIVuIIlIl^/xRS) < ||Vti(t*)|| L 2 (R 3) + ||u|| L io/3||u|| L ioJ|Vu|| L io/3 + || U || 1 ioJ| VU || L io/3 • 

Analogously, we obtain 

!MIl?LJ(/xR 3 ) ^ ll«(OIU 2 (R 3 ) + ||«||l,10.||«||^10/3 + |M|£loJM| L 10/3. 

Adding the last two inequalities and taking the supremum over all admissible pairs 
(q, r) we infer 

IMUi(/xR 3 ) < \\Vu(t*)\\ L 2 (R3) + ||«|| I ,ioJ|V'|||i + ||u|| 4 ioJ|u|| 5 i. 

Since we already know that ||u||x,io C(\\u \\ L 2, E K (0)), we can divide the time- 
interval [0, T] for any T > into subintervals Ij, j = 1, . . . , N, such that 
IHIi^^xP) < £, V j = 1,. ..,N. 

Note that N only depends on e > and on the constant C(||uo||l 2 ), E k (0)). Hence, 
for e sufficiently small a standard bootstrap argument yields 

IMU^xR 3 ) < C , (||u ||L2, J B K (0))||Vu(t i )l|L 2 (R 3 ), 
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where tj € Ij. Since we also know that ||Vtt(t)||i2 ^ C(\\uo\\ L 2, E K (0j) for each 
t > 0, we conclude 

|M| <?!(/, XR3) < C(\\U \\ L 2,E K (0)), Vj = l,...,N. 

Thus, by summing over the N subintervals we obtain the desired result. □ 

4. Proofs for subcritical damping and the time-decay of solutions 

This Section is devoted to the remaining proofs for Corollaries 2.2 - 2.5. To 
this end, we first note that in the case of an energy-subcritical damping term, 
one easily concludes local in-time well-posedness in £ by classical arguments, see 
[12, 35, 7, 8, 38]. Indeed we have the following blow-up alternative: 

Lemma 4.1. Let V be a quadratic confinement of the form (1.2) and u € S. 

Moreover, let p ^ 3, if d — 1,2 and 3 ^ p < 5, if d = 3. TTien i/iere exists a unique 
local in-time solution u G C([0,T),£). Moreover, if T < +oo, t/ien 

Urn |K*,-)||e = oo, 

In order to continue this local-in-time solution for all times we again need to 
derive suitable a-priori estimates. This will be done in the following lemma. 

Lemma 4.2. Let V be a quadratic confinement of the form (1.2) and u(t) G £ be 
a solution to (1.1). T/ien, if either: 

• p > 3, for d = 1,2, respectively 3 < p < 5, /or d = 3, 

• or p = 3 and cr ^ max {0, —A}, 

it holds 

||u(i,-)|| E <<7(||uo||j:), VOO. 

Combining this uniform bound on the energy norm of u with the assertion of 
Lemma 4.1, consequently proves Corollary 2.2 and 2.3. 

Proof. We first consider the case where p > 3 and consider the following energy- type 
functional 

E K , P (t) ■= I hvu(t,x)\ 2 + V(x)\u(t,x)\ 2 + ^u^x)] 4 + Km,x)\P +1 dx, 

JR d * 1 

with k > to be chosen later on. The time-derivative of E K>p (t) is then found to 
be 

^-E K . p (t)=-a(p-l)[ H^IV^I 2 dz-o- / luf-^Vufdx 

tit JM d JR d 

+ 7(p+l)(p- 1) / \u\ p - 3 Vp-Jdx-2a [ V\u\ p+1 dx 
-<7/c(p+l) / \u\ 2p dx-2\a [ \u\ p+3 dx. 

JR d JR d 

We rewrite, similarly as before, 

/* \u\ p - 3 W p ■ J dx = - f \u\ p ~ 3 -Vp-J dx+ / M^lViipdx, 
and also use the interpolation estimate 

(p-3)(p+l) 4p 

IMl£?. < NUT 1 IHI^ 1 . 

Using this, and following the arguments given in the proof of Proposition 3.2, we 
can always find a k = n(d,p) > such that 

E K , p {t) < E K , P (0) + c P (|KM, Vt^ 0. 
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This consequently implies a uniform bound on \\u(t, by a interpolation argu- 
ments similar to those given in the proof of Corollary 3.4. 

In a second step, we turn to the threshold situation p = 3: In the defocusing 
case A > one can use the energy-type functional specified above with k = and 
p = 3. One analogously proves that E ^(t) is decreasing along solution u{t) and 
consequently concludes that \\u(t, is uniformly bounded. 

In the focusing case A < we have to argue slightly differently. Recall that in 
this case, the NLS type equation can be written as 

(4.1) id t u = -;jA?i + V{x)u - (|A| + ia)\u\ 2 u. 

We consequently consider the corresponding linear energy functional 

Eiin(t):= [ l\Vu(t,x)\ 2 + V{x)\u{t,x)\ 2 dx. 

Differentiating Eu D (t) w.r.t time yields 

^-Ei in (t)=-[ Re(Au(t,x)d t u(t,x))dx- 2a [ V(x)\u\ 4 dx 
dt J R d J R d 

and using equation (4.1) we obtain 

E hn (t) = a [ \u\ 2 Rc(uAu) dx + |A| f \u\ 2 lm(uAu) dx - 2a f V{x)\u\ A dx 
dt J-Rd J R d J R d 



This can be re-written as 
d_ 
~dt 



-77/ \^p\ 2 dx-a \TiWu\ 2 dx + 2\X\ Re(wVu) • Im(?lVu) dx. 
ritten as 

Ey ln (t) ^ - °- [ |Vp| 2 dx - (a - |A|) f p\Wu\ 2 dx 



1/ 



^Re(u\7u) — Im(uVti) 



dx, 



and thus, if a ^ |A|, we consequently obtain E\i n (t) ^ Eu n (0) < oo. This yields 
global well-posedness of the considered NLS (4.1) and we are done. □ 

Remark 4.3. Note that the blow-up alternative given in Lemma 4.1 is not exactly 
the same as in the case of the usual NLS (without damping), for which it is enough 
to control ||Vw(£, -)||l 2 - To see this, consider the following focussing NLS (cubic, 
for simplicity, but a generalization to other power-type nonlinearities is straightfor- 
ward) : 

id t u = —-Au + V(x)u — \u\ 2 u. 

It is well known that this equation admits local in-time solutions which in addition 
preserve the mass M(t) = \\u(t)\\ L 2 and the energy functional 

E{t)= I l\Vu(t,x)\ 2 + V(x)\u(t,x)\ 2 ~hu\ 4 dx. 

JR d ^ I 

Using the Gagliardo-Nirenberg inequality 

\\u\\i^C\\u\\^ d \\Wu\\ d L ,, 

one observes that, in view of mass conservation, the (conserved) energy is a sum of 
three terms, two of which are bounded, provided that ||Vu(i, -)IIl 2 < °°- Thus, also 
the third term, i.e. the linear potential energy oc -)||l 2 , has to be bounded. 

This shows that unless ||Vu(f, -)||l 2 becomes unbounded, -)||s is a continuous 
function in time. In other words, even in the case of a given quadratic confinement, 
the possible blow-up of solutions can essentially be regarded as a local phenomena 
in x (analogously to the case without potential). 
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Our situation is a bit more involved, though: As we have seen, finite-time blow- 
up is only possible for A < and p — 3. In this case, however, we are only able to 
derive boundcdness of En n (t), provided that a < |A|. We therefore do not get any 
further insight on the nature of the blow-up and can not rule out the possibility 
that the H 1 (R 3 ) norm of the solution stays bounded but ||a;ti(i, -)IU 2 — * 00 (even 
though such a situation seems to be unlikely). 

Remark 4.4. A closely related observation concerns the following: If instead of 
(1.2), we would consider a repulsive potential of the form V(x) = — ^\x\ 2 , we would 
not succeed with our approach. More precisely, in the case of a repulsive potential, 
the corresponding term in the time-derivative of E KjP (t) comes with the wrong sign 
and thus we can not conclude that E KtP (t) is non-increasing. This is remarkable 
insofar as it is known that (sufficiently strong) repulsive quadratic potentials are an 
obstruction for the possible blow-up of solutions to NLS, see [11]. In our case, the 
situation is not so clear and global well-posedness for repulsive potentials remains 
an interesting open problem. 

Having set up global well-posedness of the equation 1.1, we finally turn to the 
proof of the decay of its solutions u(t, x) as t — > +oo. 

Proof of Corollary 2.5. Let u € C([0, oo); £) be the global solution to (1.1). From 
the dissipation equation (1.4) we obtain that 

/ / \u\P +1 dxdt^C(\\u \\ L 2). 

JO JM" 

Consider a sequence of time-steps (t n ) n &h tending to infinity and define 

u n (t,x) := u(t + t n ,x). 

From the space-time bound given above, we consequently infer that, as n —* oo: 

w„^0, in L p+1 ([0,oo) x R d ). 

On the other hand, from our global in-time existence theory, we know that (u n ) n& ti 
is uniformly bounded in C([0, oo);H 1 (W d )) and that (d t u n ) n e^ 1S uniformly bounded 
in C([0, oo); H~ 1 (M. d )). Thus by the Aubin-Lions Lemma we conclude that (u n ) n efi 
is relatively compact in C([0, T]; L 2 (R d )), for each < T < oo. We consequently 
have that there exists a subsequence (u ne ) , such that 

u ne ^0 inC([0,T];L 2 (K d )). 

From this, we consequently obtain (2.3), since 

IKWIU 2 < SU P ll u n*(*)IU 2 0- 

In addition, by using the smoothing properties of the Schrodinger group U (t) , see 
[15], we get that (u n ) ne ^ is uniformly bounded in L 2 ((0,T); H^(R d )). Invoking 
a compactness argument given in [30] we infer that (u n ) n ^ is pre-compact in 
i 2 ((0, T); H{ oc (R d )), which concludes the proof. □ 

Acknowledgement: The authors want to thank R. Carles for helpful discus- 
sions. 
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